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Accessing the connection between classical chaos and quantum many-body systems has been a 
long-standing experimental challenge. Here, we investigate the onset of chaos in periodically driven 
two-component Bose-Einstein condensates, whose small quantum uncertainties allow for explor¬ 
ing the phase space with high resolution. By analyzing the uncertainties of time-evolved many- 
body states, we find signatures of elliptic and hyperbolic periodic orbits generated according to the 
Poincare-Birkhoff theorem, and the formation of a chaotic region at increasing driving strengths. 

The employed fluctuation analysis allows for probing the phase-space structure by use of only short- 
time quantum dynamics. 
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The Poincare-Birkhoff scenario is generally recognized 
as a key mechanism for the transition from integrability 
to chaos in Hamiltonian systems HU]. It makes state¬ 
ments about the emergence of stable and unstable peri¬ 
odic orbits in the phase space of non-integrable Hamil¬ 
tonian systems. For a driven anharmonic oscillator, it 
implies that a resonant periodic orbit of the undriven 
system survives only for a pair (or more generally sev¬ 
eral pairs) of phase relations between the system and the 
driving. One of these periodic orbits is stable (more pre¬ 
cisely elliptic), which implies that a small initial deviation 
will result in quasi-periodic regular motion in its neigh¬ 
borhood. The other one is unstable (hyperbolic), which 
leads to an exponential growth of small deviations from 
this periodic orbit. This gives rise to chaotic dynamics in 
the driven system. At stronger driving, more and more 
near-resonant orbits will be affected by this resonance 
and turn into quasi-periodic or chaotic orbits, eventu¬ 
ally leading to an appreciably large chaotic region in the 
phase space of the oscillator. The controlled generation 
of such mixed phase spaces is a pathway to the creation 
of many-particle entangled states at the crossover re¬ 
gion between chaotic sea and stable islands EHS], as well 
as chaos-assisted tunneling between stable islands [Ms]. 
The experimental observation of the Poincare-Birkhoff 
scenario in many-body systems has so far been obscured 
by the large extension of the employed quantum states 
with respect to the phase-space structure |9HI3]. 

In this letter, we employ time-periodically driven Bose- 
Einstein condensates [HUS], whose small quantum un¬ 
certainties allow for exploring the phase space with high 
resolution [I6]|I7]. In this many-particle quantum sys¬ 
tem, we demonstrate signatures of elliptic and hyper¬ 
bolic periodic orbits generated according to the Poincare- 
Birkhoff theorem, and the formation of a chaotic region 


at increasing driving strengths. The employed analysis 
of the uncertainties of time-evolved many-body states al¬ 
lows us to probe the corresponding phase-space structure 
by use of only short-time quantum dynamics. 

Our system of an interacting two-component Bose- 
Einstein condensate (BEG) with linear coupling (see 
Fig. la) constitutes an anharmonic oscillator if the 
strengths of coupling Ut and interaction-induced nonlin¬ 
earity Nx are comparable. The phase and amplitude 
controlled coupling in our system is implemented using a 
two-photon transition between the two hyperfine states 
i;) =\F = l,mF = 1) and |t) =| 2, -1) of ^^Rb, and the 
nonlinearity is enhanced in the vicinity of a Feshbach res¬ 
onance (details in the Supplemental Material [H]). Our 
experiments investigate the phase space spanned by the 
population imbalance ^ = {N^ — Ni)/N and the relative 
phase ip = {cpI — cp^) between the two states. Here, 
and are the atom numbers of the individual compo¬ 
nents, which constitute the experimental observables (see 
Fig. la), and N « 700 is the total atom number. External 
dynamics, central to related proposals and experiments 
with driven BECs [HI [TSl US] , are frozen out and the 
dynamics happens exclusively in the collective internal 
degree of freedom. The extension of a minimal quantum 
uncertainty state in this system, i.e. a coherent spin state 
(CSS), is given by [z, cp] = PljN, which can be interpreted 
as realizing an effective Planck’s constant /leff = 2/A/' (see 
Fig. lb). While the corresponding uncertainty is large for 
few-state systems m, /ieff in our mesoscopic two-mode 
BEC is small and allows for probing the phase space with 
high resolution. 

For N ^ oo, the dynamics of our system can be de¬ 
scribed by the classical Hamiltonian [4] 

H{z, p,t) = ^z‘^ — \/l — z‘^ cos p P ez (1) 

2 Hq 
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FIG. 1. (Color online). Poincare-Birkhoff scenario in a many-body quantum system. (a) A two-component Bose-Einstein 
condensate (N ^ 700 atoms) with interspecies interaction x and linear coupling Q realizes an anharmonic oscillator with the 
population imbalance and the relative phase (f as the canonical variables. State selective absorption imaging allows for the 
precise determination of the populations of the two components, (b) The minimal uncertainty state is a coherent spin state 
whose extension in the phase-space representation (z,(p) corresponds to hefi = 2/N. (c) Classical orbits for the undriven case 
{A = 0) at A = Nx/^o = 0.7 (green square: elliptic fixed point at (f = tt). (d) Stroboscopic Poincare sections, showing 
the formation of an alternating sequence of elliptic (blue triangles) and hyperbolic (red dots) periodic orbits according to the 
Poincare-Birkhoff theorem in the case of weak driving (A = 0.03), the emergence of a chaotic layer for A = 0.2, and a larger 
chaotic region at even higher modulation amplitude. 


where the nonlinearity A = Nx/^o and the detuning 
e = are rescaled with the Rabi frequency Qq [18]. 

The drive is realized by 0(t) = Qq {1 A sin [uj{t to)]). 
The resulting phase space is represented on a sphere with 
the polar angle = arccos 2 : and the azimuthal angle 
taking into account that 2 : and ip are bounded. In Fig. Id, 
it is visualized for different driving amplitudes by stro¬ 
boscopic Poincare sections at times which are multiples 
of T = The undriven system for A < 1 features 
an elliptic fixed point near [z = 0,(p = n) (Fig. Ic) and 
anharmonic oscillations around it. By driving with a 
frequency cj, a 2 : 1 resonance is realized for the orbit 
with the unperturbed frequency uj/2, corresponding to 
a specific amplitude of the anharmonic oscillator (exam¬ 
ple for uj = l.bflo shown in Fig. Id). According to the 
Poincare-Birkhoff theorem, this implies the destruction 
of the resonant orbit except for two elliptic and two hy¬ 
perbolic periodic orbits. The resulting stationary points 
in the Poincare section are marked as blue triangles and 
red stars. For increasing driving, the phase space acquires 
more and more structure, leading to a mixed phase space 
characterized by regular islands and a chaotic sea (see 
Fig. Id). Similar behavior has also been theoretically in¬ 
vestigated for spin-exchange dynamics in Bose-Einstein 
condensates [25ll26] . 

For experimentally revealing this Poincare-Birkhoff 
scenario, we prepare an initial coherent state close to 
either an elliptic or a hyperbolic periodic orbit and ob¬ 
serve the subsequent time evolution of mean and variance 
of the population imbalance z. During the time evolu¬ 
tion, the elliptic and hyperbolic periodic orbits perform 
a counterclockwise rotation around the central elliptic 
fixed point, as shown in Fig. 2a. This implies identical 
phase spaces after one driving period T, and a full rota¬ 
tion of elliptic and hyperbolic periodic orbits after 2T. 


All phase-space structures, e.g. regular orbits close to the 
elliptic periodic orbits as well as the chaotic layer, follow 
this rotation. 

We initialize our system with a CSS centered at {zq = 
0.55, (/:?o = tt), which is located close to the hyperbolic 
periodic orbit for to = 0.2 T (red) and close to the ellip¬ 
tic periodic orbit for to = 0.6 T (blue) for the parameters 
A = 0.2, A = 0.7, e = -0.11 and u; = l.bllo* The 
mean of the quantum state (Fig. 2b) oscillates in con¬ 
junction with the hyperbolic (red) or the elliptic orbit 
(blue), and resembles the dynamics of the undriven sys¬ 
tem. However, the measured variances, corresponding to 
the 2 : projection of the quantum uncertainty of the time- 
evolved states, show striking differences, as depicted in 
Fig. 2c. In classical dynamical systems, two initially adja¬ 
cent points at the unstable manifold or within a chaotic 
region are expected to separate exponentially in time. 
Thus, for the preparation close to the hyperbolic orbit, 
the normalized variance [18] exhibits an overall increase 
compared to the undriven case. The observed additional 
oscillations result from the spreading of the distribution 
along its unstable manifold, which rotates around the 
central elliptic fixed point of the undriven system (inset of 
Fig. 2c). Thus it constantly changes the orientation to the 
readout direction This behavior is well reproduced by 
propagating the classical Hamiltonian with an ensemble 
of initial conditions corresponding to a CSS (black solid 
lines in Fig. 2, [18]; quantum calculations for N = 700 
show no significant deviations). The numerical calcula¬ 
tions underestimate the fluctuations, since experimental 
imperfections such as magnetic field fluctuations, parti¬ 
cle loss and the associated decoherence and change of the 
system parameters are not included. 

The initial CSS close to the elliptic periodic orbit re¬ 
mains bounded and even shows a smaller variance com- 
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FIG. 2. (Color online). Detecting the Poincare-Birkhoff scenario via variance analysis, (a) Stroboscopic Poincare sections for 
A = 0.7 and A = 0.2 for various times t G [0,2T], showing one full counterclockwise rotation of the elliptic and hyperbolic 
periodic orbits around the central elliptic fixed point. Red and blue points: Time evolution for {zq = 0.55, (/?o = tv) and driving 
time offsets to = 0.2T (red) and to = 0.6T (blue), (b) The experimentally obtained mean values of exhibit oscillations with 
period 2T for preparations close to the hyperbolic orbit (left, red) and the elliptic orbit (right, blue). Error bars are the RMS 
deviation. Black lines and shaded regions: Mean and RMS uncertainty for classical time evolution of an ensemble corresponding 
to a CSS. (c) Variances of in units of the corresponding coherent state m- Preparation near the hyperbolic orbit (left) results 
in a strong increase compared to the undriven dynamics (gray), whereas suppressed uncertainties can be observed close to the 
elliptic orbit (right, blue). The results of classical simulations are shown as black lines. The oscillations of the variance result 
from the rotation of the corresponding manifolds around the central elliptic fixed point (insets). 


pared to the undriven case (Fig. 2c right). This obser¬ 
vation builds on the size of a Planck cell in our system, 
which is small compared to the phase-space region with 
quasi-periodic oscillation about the elliptic periodic or¬ 
bit. This is the many-particle analogon of non-dispersive 
single-particle wave packets, which were predicted and 
observed in microwave driven Rydberg atoms [27H29] . 
Such a wave packet retraces the periodic motion of the 
elliptic periodic orbit without being subjected to spread¬ 
ing. 

We can access a mixed phase space with regular and 
chaotic regions by changing the relative nonlinearity to 
A = 1.5, the detuning to e = —0.07 and the driving 
strength to A = 0.07. For this scenario, the phase space 
in the absence of driving exhibits one hyperbolic and two 
elliptic fixed points (Fig. 3a, top). The driven system 
shows, embedded in a large chaotic region, two regu¬ 
lar islands originating from the two elliptic fixed points 
of the undriven system and two moving regular islands 
emerging due to a 2:1 resonance for a driving frequency 


uj = 1.6 flo- Here, we focus on the differences in the 
time evolution of states that are initially prepared close 
to this moving island or within the chaotic region. For 
fixed initial preparation {zq = 0, (^o = 2.51), we vary 
the time offset to in order to start the dynamics within 
the chaotic region {to = 0.9T) or close to the moving 
island {to = OAT). We access two axes of the phase 
space by measuring particle number difference with and 
without a final 7r/2 pulse, corresponding to readout of 
the 2 : and the cp direction. The undriven system shows a 
moderate broadening of the distribution along both axes 
(Fig. 3, upper panels). In the driven system, starting 
in the chaotic region we find fast broadening, eventually 
covering large parts of the phase space along the z di¬ 
rection. In contrast, by starting in the vicinity of the 
moving island, the distribution remains confined within 
the observation time and disperses even less than in the 
undriven system m- This behavior can be intuitively 
grasped from the structure of the corresponding phase 
space (chaotic vs. sticky region) and is in good agreement 
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FIG. 3. (Color online). Signatures of a mixed phase space in 
the short-time quantum evolution, (a) The undriven system 
(A = 1.5) exhibits one hyperbolic and two elliptic fixed points 
at (/? = TT. For A = 0.07, the strongly driven regime with a 
mixed phase space is reached, which is probed by different 
initial preparations (depicted by 2a areas) with the time off¬ 
sets to = 0.9T (center, red) and to = 0.4T (bottom, blue), (b) 
Time evolution of the z distributions in gray scale. Compared 
to the undriven case (top), we observe increased fluctuations 
in the case of a preparation deep in the chaotic sea (mid¬ 
dle) and reduced fluctuations close to a moving island (bot¬ 
tom). Mean and RMS deviation of the classical simulations 
are shown as solid and dashed lines, (c) Similar signatures of 
the underlying phase space are observed in y = sin((/?). 


FIG. 4. (Color online). Due to the small size of heff in our 
system, we can experimentally probe phase-space structures 
on small scales. The initial position is adjusted by varying the 
time offset to of the driving, and is indicated in the Poincare 
section for to = T along with the size of the effective Planck 
constant (left, blue points with 2a area in red), (a) Weakly 
perturbed system: After a fixed evolution time (t = 48 ms), 
we hnd enhanced variances near the hyperbolic periodic orbit 
and reduced variances near the elliptic orbit. Our hndings 
are compared to the undriven case (gray) and the classical 
simulation (black line), (b) Strongly perturbed regime: The 
variance after an evolution time of t = 61 ms reveals a strong 
dependence on the initial preparation, indicating the high res¬ 
olution in phase space. 


with the predictions from classical simulations. Note that 
in the simulations the initial state is set to (zq = —0.3, 
(po = 2.68), which captures loss effects during the exper¬ 
iment [18]. 

Fig. 4 shows the results of our investigation of the 
transition from regular to chaotic regions for the weakly 
(Fig. 4a) and the strongly perturbed system (Fig. 4b). In 
both cases, we observe a clear correspondence of the fi¬ 
nal variance after a fixed time to the structure of the 
phase space at the initial preparation. In the weakly 
perturbed system, the variance increases in the vicinity 
of the hyperbolic periodic orbit and is minimal (com¬ 
pressed compared to the undriven case) in the regular 
region close to the elliptic periodic orbit. A similar be¬ 
havior is found for the strongly perturbed system with 
a strong increase of variance in the chaotic region and a 
sudden drop when approaching the sticky regions close 
to the moving island. This drop occurs within 3.5 widths 
of a minimal uncertainty state (CSS) and demonstrates 


the high experimental resolution of the phase space. The 
observations are in good agreement with the simulation 
(Fig. 4, solid black line). 

Our experiments reveal that the mixed phase-space 
structure of a driven system can be probed even in the 
short time dynamics of just a few driving periods by 
a variance analysis of the quantum states. Further¬ 
more, the experimental generation and verification of the 
Poincare-Birkhoff scenario is a cornerstone to control the 
dynamics of complex many-particle systems. Since the 
phase-space structure, /ieff and decoherence processes in 
the system can be precisely controlled, it is a model sys¬ 
tem for the exploration of the classical limit of quantum 
many-body systems 13 EDI ED. 

We are grateful to S. Gardiner and C. Weiss for 
discussions and valuable feedback on the manuscript. 
We thank P. Kevrekidis and L. Carr for revealing con¬ 
nections to general nonlinear dynamics and D. Linne- 































5 


mann for experimental assistance. We acknowledge fi¬ 
nancial support through the DFG Forschergruppe 760 
’Scattering Systems with Complex Dynamics’ of the 
Deutsche Forschungsgemeinschaft, the Heidelberg Cen¬ 
ter for Quantum Dynamics and the European Comis¬ 
sion FET-Proactive grant AQuS (Project No. 640800). 
W.M. acknowledges support by the Studienstiftung des 
deutschen Volkes. 


PoincareBirkhoff@matterwave.de 

[1] G. D. Birkhoff, On the periodic motions of dynamical sys¬ 
tems. Acta Mathematica 50 , 359-379 (1927). 

[2] E. Ott, Chaos in dynamical systems. (Cambridge Univ. 
Press, Cambridge, 2002). 

[3] C. Weiss and N. Teichmann, Differences between Mean- 
Field Dynamics and N-Particle Quantum Dynamics as a 
Signature of Entanglement. Phys. Rev. Lett. 100 , 140408 
(2008). 

[4] Q. Xie and W. Hai, Quantum entanglement and chaos 
in kicked two-component Bose-Einstein condensates. Eur. 
Phys. J. D 33, 265-272 (2005). 

[5] S. Chose, R. Stock, P. lessen, R. Lai, and A. Silber- 
farb. Chaos, entanglement, and decoherence in the quan¬ 
tum kicked top. Phys. Rev. A 78 , 042318 (2008). 

[6] W. A. Lin and L. E. Ballentine, Quantum Tunneling and 
Chaos in a Driven Anharmonic Oscillator. Phys. Rev. 
Lett. 65, 2927-2930 (1990). 

[7] S. Tomsovic and D. Ullmo, Chaos-assisted tunneling. 

Phys. Rev. E 50 , 145-162 (1994). 

[8] M. P. Strzys, E. M. Graefe, and H. J. Korsch, Kicked 
Bose-Huhhard systems and kicked tops - destruction and 
stimulation of tunneling. New J. Phys. 10 , 013024 (2008). 

[9] F. L. Moore, J. C. Robinson, C. Bharucha, P. E. Williams, 
and M. G. Raizen, Observation of Dynamical Localization 
in Atomic Momentum Transfer: A New Testing Cround 
for Quantum Chaos. Phys. Rev. Lett. 73 , 2974-2977 
(1994). 

[10] F. L. Moore, J. C. Robinson, C. F. Bharucha, B. Sun- 
daram, and M. G. Raizen, Atom Optics Realization of the 
Quantum 6-Kicked Rotor. Phys. Rev. Lett. 75 , 4598-4601 
(1995). 

[11] D. A. Steck, W. H. Oskay, and M. G. Raizen, Observation 
of Chaos-Assisted Tunneling Between Islands of Stability. 
Science 293 , 274-278 (2001). 

[12] W. K. Hensinger et al.. Dynamical tunnelling of ultracold 
atoms. Nature 412 , 52-55 (2001). 

[13] S. Chaudhury, A. Smith, B. E. Anderson, S. Ghose, and 
P. S. lessen. Quantum signatures of chaos in a kicked top. 
Nature 461 , 768-771 (2009). 

[14] G. 1. Duffy, S. Parkins, T. Muller, M. Sadgrove, R. Leon¬ 
hard!, and A. C. Wilson, Experimental investigation of 
early-time diffusion in the quantum kicked rotor using 
a Bose-Einstein condensate. Phys. Rev. E 70 , 056206 
(2004). 

[15] G. Behinaein, V. Ramareddy, P. Ahmadi, and 
G. S. Summy, Exploring the Phase Space of the Quan¬ 
tum 5-Kicked accelerator. Phys. Rev. Lett. 97 , 244101 
(2006). 

[16] H. Frahm and H. 1. Mikeska, On the Dynamics of a 


Quantum System which is Classically Chaotic. Z. Phys. B 
Condensed Matter 60 , 117-126 (1985). 

[17] G. L. Salmond, C. A. Holmes, and G. 1. Milburn, Dy¬ 
namics of a strongly driven two-component Bose-Einstein 
condensate. Phys. Rev. A 65, 033623 (2002). 

[18] See Supplemental Material for details on the experimen¬ 
tal sequence, data analysis and simulations. 

[19] A. Micheli, D. laksch, 1. 1. Cirac, and P. Zoller, Many- 
particle entanglement in two-component Bose-Einstein 
condensates. Phys. Rev. A 67 , 013607 (2003). 

[20] T. Zibold, E. Nicklas, C. Gross, and M. K. Oberthaler, 
Classical Bifurcation at the Transition from Rabi to 
Josephson Dynamics. Phys. Rev. Lett. 105 , 204101 
( 2010 ). 

[21] H. Strobel et al., Fisher information and entanglement 
of non-Caussian spin states. Science 345, 424-427 (2014). 

[22] S. A. Gardiner, D. Jaksch, R. Dum, J. 1. Cirac, and 
P. Zoller, Nonlinear matter wave dynamics with a chaotic 
potential. Phys. Rev. A 62 , 023612 (2000). 

[23] T. S. Monteiro, A. Rancon, and J. Ruostekoski, Nonlin¬ 
ear Resonances in 5-Kicked Bose-Einstein Condensates. 
Phys. Rev. Lett. 102 , 014102 (2009). 

[24] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, 
Quantum Coherent Atomic Tunneling between Two 
Trapped Bose-Einstein Condensates. Phys. Rev. Lett. 79 , 
4950-4953 (1997). 

[25] J. Kronjager, K. Sengstock, and K. Bongs, Chaotic dy¬ 
namics in spinor Bose-Einstein condensates. New J. Phys. 
10 , 045028 (2008). 

[26] J. Cheng, Chaotic dynamics in a periodically driven spin- 
1 condensate. Phys. Rev. A 81 , 023619 (2010). 

[27] A. Buchleitner, D. Delande, and J. Zakrzewski, Non- 
dispersive wave packets in periodically driven quantum sys¬ 
tems. Physics Reports 368 , 409-547 (2002). 

[28] H. Maeda and T. F. Gallagher, Nondispersing Wave 
Packets. Phys. Rev. Lett. 92 , 133004 (2004). 

[29] H. Maeda, D. V. L. Norum, and T. F. Gallagher, Mi¬ 
crowave Manipulation of an Atomic Electron in a Classical 
Orbit. Science 307 , 1757-1760 (2005). 

[30] S. A. Gardiner, J. 1. Cirac, and P. Zoller, Quantum Chaos 
in an Ion Trap: The Delta-Kicked Harmonic Oscillator. 
Phys. Rev. Lett. 79 , 4790-4793 (1997). 

[31] S. Habib, K. Shizume, and W. H. Zurek, Decoherence, 
Chaos, and the Correspondence Principle. Phys. Rev. 
Lett. 80 , 4361-4365 (1998). 



6 


SUPPLEMENTAL MATERIAL: OBSERVING THE EMERGENCE OF CHAOS IN A MANY-PARTICLE 

QUANTUM SYSTEM 


Experimental system 

For our experimental studies, we employ BECs con¬ 
taining N ^ 700 atoms in the two hyperfine levels 
\D =\ F = l.rriF = +1) and |t) =\ F = 2,mi? = — 1). 
In the presence of linear coupling and collisional inter¬ 
action, this system can be described with collective spin 
operators and is governed by the Lipkin-Meshkov-Glick 
Hamiltonian H = xJz ~ ^Jx F 5Jz mm- 

The coupling is realized via a two-photon microwave 
and radio frequency (RE) transition with Rabi frequency 
Q. A slight detuning 6 of the coupling from the transi¬ 
tion breaks the symmetry of the phase space along the 2 ; 
direction. The coupling phase, frequency and amplitude 
can be controlled using an arbitrary waveform generator 
for the RE, and the amplitude of the microwave radia¬ 
tion can be changed by use of an attenuator on a two-port 
switch. Eor all experiments, AC Zeeman frequency shifts 
resulting from the two-photon coupling scheme and col- 
lisional mean field shifts are compensated. 

An interspecies Eeshbach resonance at a magnetic field 
of 9.1 G permits the enhancement of the atomic nonlin¬ 
earity to Nx ~ 27r X 32(2) Hz at our working field of 
9.12 G. The resulting atom loss during the evolution time 
(typical time scale 80 ms) causes a time-dependence of 
the atom number dependent collisional shifts and thus 
of the detuning 6. This shift is compensated by an inde¬ 
pendently calibrated feed forward on the frequency of the 
RE radiation. Additional temporal changes of x(A^) are 
not compensated, but do not change the overall structure 
of the phase space. Technical fluctuations and drifts in 
the detuning resulting from magnetic held changes are as 
small as as ~ 27r x 0.45 Hz, as deduced from independent 
Ramsey measurements. 

An optical lattice potential is used to simultaneously 
generate about 30 independent condensates, which yields 
enhanced statistics and freezes out external degrees of 
freedom due to the high trapping potential in all spatial 
directions. The individual condensates are imaged using 
site- and state-resolved absorption imaging after Stern- 
Gerlach separation of the two components. 

The classical limit of this system is obtained for N 
00 and keeping A constant by replacing the operators 
with their mean values. The resulting mean held dy¬ 
namics have been intensely studied both theoretically 
and experimentally mm- In this classical description, 
a bifurcation [2] occurs for similar strength of coupling 
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and nonlinearity at Acr = ^l + |e|i^^, where one ellip¬ 
tic fixed point turns into one hyperbolic and two elliptic 
fixed points. 


Experimental sequence 

At the beginning of the experimental sequence, all 
atoms are initialized in |^). A high-power coupling pulse 
(O « 27r X 310 Hz) is employed to prepare the initial quan¬ 
tum state at a fixed initial population imbalance zq by 
properly adjusting the length of the pulse. 

Subsequently, the power of microwave and RE radia¬ 
tion is attenuated to reach the desired value of A, and the 
phase of the RE radiation is nonadiabatically adjusted to 
prepare the initial state at the chosen phase cpo. The driv¬ 
ing is realized via the corresponding periodic modulation 
of the RE signal amplitude. After the evolution time t, 
the two components are absorption imaged to obtain the 
^ projection of the final state. Alternatively, for readout 
along the y = sin((/9) direction, we apply a further fast 
7r/2 rotation around the x direction, which converts the 
phase information into population imbalance. 

Data analysis 

The results are postselected on atom number following 
the actual loss rate in order to ensure a consistent en¬ 
semble starting with N = 700 ± 25 at t = 0. Each data 
point comprises 60 to 90 independent realizations of the 
experiment. Outliers are rejected by a modified Z-score 
test. 

The obtained variances in z are normalized by 
Varcss = (1 “ (^)^) x where N{t) denotes the to¬ 

tal atom number at time t. This ensures that the state 
is referenced to the corresponding CSS regardless of its 
position on the sphere or evolution time in the presence 
of atom loss. The error bar of the variance is given by 
the 1 s.d. confidence interval obtained from Jackknife 
resampling. 

Simulations 

Eor the comparison of the many-particle quantum 
system with the classical model (Eq. 1 of the main text), 
the finite width of the wave function is accounted for by 
integrating the corresponding equations of motion for 
a set of points sampled from the Gaussian distribution 
of the corresponding coherent spin state. Erom these 
integrations, the corresponding mean values and vari¬ 
ances are obtained. While the change of the collisional 
shifts due to atom loss is experimentally compensated, 
the corresponding change of the relative nonlinearity 
A is not corrected. To obtain a proper phase space 
description, this time dependence is captured in the 
simulations by using an effective starting point, which is 
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obtained from fits for both parameter sets from reference 
measurements with the undriven system. 
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